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Abstract 

A 125 GeV Higgs renders the simpler GMSB models unnatural, essentially pushing the soft spectrum beyond the 
LHC reach. A direct coupling of the matter and messenger fields, that facilitates an enhanced mixing in the squark sector, 
is a way to ameliorate this deficiency. We construct all possible messenger-matter interaction terms considering the 
messenger multiplets in 1, 5 and 10 dimensional representations of the SU(5). A Froggatt-Nielsen like flavor framework 
connected with the origin of fermion mass hierarchy is utilized to control the interaction terms and suppress FCNC. We 
perform a detailed comparative study of the efficiency of such interaction terms to boost the Higgs mass keeping the soft 
spectrum light. We identify the more promising models and comment on their status in present and future collider studies. 



1 Introduction: 
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Identification of the recently observed 0]|2] scalar field at the LHC with the Higgs would cast a long shadow on Gauge 
Mediated Supersymmetry Breaking (GMSB) models Q0BI- This is a direct consequence of the well known phenomenon 
that in pure gauge mediation models, mixing in the scalar sector is minimum. This implies that in order to raise the Higgs 
mass to ~ 125 GeV one needs stop masses at several TeV, castigating these models to an unnatural existence with bleak 
chances to be probed at collider experiments like the LHC. 

A resolution of this predicament is to consider direct coupling between the messenger and matter fields or models with 
gauge messengers 0. These scenarios can in principle generate a sizable trilinear coupling that can boost the Higgs mass 
given by, 
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where, X t = A t — \i cot (3 and M s = ^m^m^ , while keeping the scalar spectrum within the range of interest for 
collider physics. In case of gauge messengers the trilinear couplings are proportional to gauge charges. However real- 
istic models recently studied in (6) predict a relatively heavy spectrum. In this paper we will consider models of direct 
messenger-matter interactions that can lead to relatively large trilinear coupling and a considerably light soft spectrum 
making them more pleasing from a fine-tuning point of view and more interesting phenomenologically. These interac- 
tion terms generally lead to new contributions to the scalar masses, thus producing a correlated perturbation of the pure 
gauge mediation soft spectrum. The strength of the interaction and hence the size of the trilinear coupling is principally 
constrained from the considerations to keep the scalar masses non-tachyonic and achieve radiative electroweak symmetry 
breaking (REWSB). Also a cause for concern are the complications in the flavor sector 0^9). These tend to push the 
messenger scale upward into regions where the Gravitino mass goes beyond ~ sub-KeV range putting unfavorable upper 
bounds on the reheating temperature inviting strong constraints from BBN iflOl . The flavor constraints are severe for the 
first two generations of fermions and can be minimized by considering that the messengers preferentially couple to the 
third generation, which is relevant for enhancing the Higgs mass. This can be ensured by imposing judiciously chosen 
flavor symmetries. 
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Several models of GMSB augmented with messenger-matter interactions to alleviate the problem with the heavy Higgs 
have been suggested in the literature lfTTHT6l . The basic idea being the coupling of the messengers to one of the Q, U, H u 
MSSM multiplets, generating a sizable A t at the messenger scale. Possible interaction terms in a given scenario get 
determined by the content and quantum numbers of the messenger sector. In this paper we make a systematic study 
of messenger-matter interactions, possible within a well defined framework. As an organizing strategy we consider the 
messengers as vector pairs embedded in one of the simpler representations of the GUT group SU(5). We will restrict 
ourselves to messengers in the 1, 5 and 10 dimensional representations of SU(5) or an admixture. This ensures that 
the perturbative unification of gauge coupling is not ruined due to introduction of these additional fields Q. For the first 
time we make a quantitative comparison of various models, including some that have been discussed in the literature and 
some entirely new scenarios, in term of their effectiveness to raise the Higgs mass without the usual pitfall of large scalar 
masses. In order to make the comparison, we numerically scan all parameters of every given model over a suitable range 
and project the allowed regions on a common parameter space. This allows us to make precise statements about scenarios 
preferred by the recent data on the supposed Higgs mass. 

For simplicity, in our study we will not attempt to model the hidden sector and simply assume supersymmetry is broken 
by the vev of a spurion field which couples to the messengers. The crucial point would be that the messengers other than 
having usual gauge couplings to the MSSM sector, now also couple directly through the superpotential. In this paper, 
we enumerate possible messenger-matter interaction terms allowed by a given messenger sector. For each scenario we 
compute the contributions of these new terms to the scalar masses at one and two loop order. We find that the one loop 
contributions are always tachyonic but they are suppressed by x 2 = (F/M 2 ) 2 where M and F are messenger scale and 
supersymmetry breaking scale respectively. Thus for a given A = F/AI or soft scalar mass, this contribution decouples 
as the scale of supersymmetry breaking is increased. The sign of the two loop contribution is model dependent, however 
there is no suppression from the supersymmetry breaking scale. In most regions of the parameter space this becomes 
the dominant contribution from the new terms. One loop renormalization group equations are used to run these soft 
parameters down from the scale of supersymmetry breaking to the weak scale and the sparticle spectrum is generated. 

The rest of the paper is organized as follows. In Section 2 we describe our organizing principle and enumerate the 
possible interaction terms. In Section 3 we compute the soft masses and trilinear couplings for the different models. In 
Section 4 we describe our numerical procedure and present our results comparing the models. Finally we conclude. 



2 Messenger-Matter interactions 

In this section, we collect the possible interactions between the messenger sector and the MSSM fields. A useful way to 
organize these is to consider the fields embedded in a representation of some GUT group like SU(5). As usual, the MSSM 
fields can be embedded into representations of the SU(5) group as follows, 

5 = (3, 1, -) ffi (1, 2, ~) = D c © L, 10 = (3, 2, -) © (3, 1, -\) © (1, 1, 1) = Q © U c © E c 
6 2 b 6 

5h = (integrated out fields) © Ha, 5h = (integrated out fields) © H u 

The messenger fields are in a vector like pair embedded in 1, 5, 10 and their conjugate representations of SU(5). For 
the rest of this paper we use the following nomenclature for the messenger sector, 

l?n = S rn , 5 TO = D m ffi H u , 5 rn = © 

10 m = Qm ffi t^ m © -E^i 10 m = Q m ffi U m ffi E^, 

where the subscript m identifies the messenger fields. 

Technically the singlet is not a valid gauge messenger in the usual sense. However, it can couple directly to the visible 
sector through superpotential terms and thus will be considered in the following discussion. Within this framework for 

'Note that it is possible to ensure perturbative gauge coupling unification without considering complete representations of the GUT groups 1171 . 
These magic combinations can potentially lead to interesting phenomenological scenarios 1181 , we will not discuss them in this paper. 
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Table 1 : U( 1) flavor symmetry and the R-parity charges of different MSSM multiplets. The subscripts denote the flavor index. Values 
p — 0, 1,2 approximately explain the fermion mass hierarchy. The charge assignment for the MSSM multiplets will remain same 
throughout the paper. 

the messenger and the matter sectors, it is simple to write down all the possible SU(5) invariants that can be constructed: 

Only Singlets (i) 5 ff 5 ff l TO (ii) 5 H 5 1 TO 

Only 5 © 5 (i) 10 5 H 5 ro (ii) 10 5 5 m (Hi) 10 5 ro 5 ro (to) 10 5 mi 5 mi (w) 10 10 5 m 
Only 10 ©TO (i) 1010 m 5# (ii)10 m 10 m 5 H (Hi) T0 m T0 m 5 H (w)T0 TO T0 m 5. 

Note that the invariants 10 m 10 m 5ff and 10 m 10 m 5 will not give A t hence will not be considered further. The possibility 
that the messenger sector is composed of fields in different representation can also exist. In this case we expect the usual 
invariants listed above should reappear. The new possible invariants are given below, 

Singlet+5©5 (i) 5#5 m l m (ii) 5 m 51 m (in) 5 m 5#l m 
Singlet + 10 © TO (i) 10101 m 
5© 5 + 10 ©TO (i) 1010 m 5 m (ii)TO m 5 m 5 H 

Again we find that terms like 5 m 5S and 5 m 5# .S will not give A t . To keep the discussion tractable we will utilize prudently 
chosen flavor symmetries to suppress interaction terms other than the ones listed above. 

As emphasized earlier a host of issues with flavor including FCNC can be controlled if we consider scenarios where 
the messenger sector preferentially couples to the third generation multiplet. An economical proposition is to connect this 
to the flavor symmetry that is responsible for the hierarchy in measured mass of the Standard Model fermions iflZl . For 
example one can consider the Froggatt-Nielsen |[T9l framework that necessitates an U(l)p flavor symmetry group, under 
which the MSSM chiral multiplets may be charged. This flavor symmetry is spontaneously broken at some high scale 
Mstring by the vev of the flavon field (Z), with a conventionally assigned flavor charge — 1, For this choice, operators 
with a negative U(1)f charge mismatch (AF) are considered absent in the effective theory. However, if AF > 
for any superpotential term then it is suppressed by the usual Froggatt-Nielsen factor e AF where e = (Z) /M s t r ing- 
Flavor symmetries of this kind can successfully explain observed fermion mass matrices within present experimental 
uncertainties. Interestingly the anomalous nature of the symmetry can motivate compensatory exotic particles that can be 
probed at present and future collider experiments ll20H22l . thus providing an handle to probe the existence of these flavor 
structures. Together with R-parity assignments this will be enough to determine the interaction terms uniquely for most 
of the scenarios that will be explored in this paper. The usual charges for the MSSM and the flavon multiplets are given 
in Table Q] 

Another symmetry that is usually useful in safely segregating the messenger and visible sectors in usual GMSB models 
is the messenger parity 0231 under which all the messenger multiplets are assumed odd while the MSSM multiplets are 
considered even. One can classify the interaction terms introduced above into two categories with different implications 
for possible messenger parities. In models where the interaction term involves two messenger multiplets and one MSSM 
multiplet one can impose messenger parity consistently. This immediately forbids mass mixing terms between the mes- 
senger and matter multiplets at all orders of the perturbation theory. While for models where the interaction term involves 
one messenger multiplet and two MSSM multiplets, a consistent messenger parity cannot be constructed. This will lead to 
messenger-matter mixing either at tree level or at higher orders, leading to non-trivial contributions to the soft spectrum. 
We will discuss the consequences of this for specific scenarios in the next section. 
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3 The soft breaking masses 



We assume that supersymmetry is broken due to some hidden sector dynamics that can be parametrized into a spurion (X) 
vev. The messengers which are charges under the MSSM gauge groups, couple to these spurion fields and communicate 
supersymmetry breaking to the visible sector through the usual gauge couplings. Usually the messengers are assumed to 
have some messenger parity that prevents these fields from mixing or interacting with the visible sector. In this paper we 
will relax this and consider all possible interaction terms between the visible and hidden sectors. The entire superpotential 
can be schematically written as, 



w = w M ssm(4>vis) + ^X4> m 4> m + w mt (fa n ,fa ls ), 



(2) 



where the spurion gets a supersymmetry breaking vev (X) = M +8 2 F and {4> m , <f> m } are vector like pair of the messenger 
fields and fa,i s are the usual MSSM chiral supermultiplets. We define the messenger scale M mess = X^M and A = F/M. 
The loop integrals can be expressed in terms of the dimensionless parameter x = A/M mess . 

The most general messenger sector however can involve mass terms for the messenger fields of the form rriijfafa . This 
can introduce a new dimensionful parameter other than the messenger scale. The ensuing complication in determination 
of the soft spectrum through the wave-function renormalization technique lead to the so called m# u — A u problem lfl3l 
in models where H u and/or Hd couple with messengers. In the models discussed in this paper we will find the these 
mass terms are either absent or suppressed by a factor e a , where a is the flavor charge of the spurion field X. We can in 
principle make a large to suppress these terms effectively. We will neglect these terms in our calculations below. 

The soft masses get the usual contributions from the gauge interactions of the messenger fields at one loop for the 
gaugino masses and two loops for the scalar masses, given by H, 
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where d is the Dynkin index and is 1 for messengers in the 5 © 5 and 3 for messengers in 10 © 10. The C r (i)'s are the 
usual Casimir invariants for the representation i and, 



/(*) = 



l + x 



log(l + x) - 2Li 2 (x/[l + x]) + -Li 2 (2x/[1 + a;]) 



-x), 



g(x) = — [(l + x)log(l + x) + (l-x)log(l-x)] 



(4) 



Now we turn to the contribution of the interaction term in Eq. [2] In this paper we will only consider interaction terms 
of the form Wi n t = Xijkfafa<f>k> where at least one chiral multiplet from both the messenger and the visible sector are 
present. The contribution at one loop level to a field fa belonging to the MSSM can be directly computed and is given by, 

^ M il- loop = -C ijfe ^pA 2 M*), (5) 

where djk is the multiplicity factor that measures the effective number of messenger fields that the MSSM supermultiplet 
fa couples to through the superpotential coupling A«fc including the appropriate group theoretic factors^. The expression 
of h(x) is as follows II13II151 , 

O (x-2)log(l-x)-(x + 2)log(l+x) 
h(x) = 3 . (6) 

a; 4 

Let us make the following observations regarding the one loop contribution: (i) note that h(x — >• 0) — s- 1 and thus (ii) 
the one loop contribution decouples with the messenger scale ( SMf (x — > 0) 1 1 _ loop —> 0), (iii) the contribution is always 
negative. 



2 For instance, in the interaction XQU^Hjj , Cijk for Q is 1 whereas Cijk for Hu is 3 (color factor). In contrast for XQmU^Hjj , Cijk for Hu 
is now 6. 
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It is easier to use the wave function renormalization techniques [24, 25) in order to compute the the two loop contribu- 
tions. They can be written in terms of the anomalous dimension and the f3 functions for the gauge coupling, below and 
above the messenger scales. Adaptations of the generic framework for messenger-matter interaction terms were computed 
in 01 111251 . we quote them for the sake of completeness, 



| 2 — loop 
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A 2 , 
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where, AX 



M„ 



A ii*\\-ioo P = ~ (7i + A 7j + Mk) h A (x) A, 
[X+ - X~]\ Mmess = [X(M mess + 5)- X(M mess - 5)]\s-*>, Px = dX/dt and 



hA(x) = l-l Ji±l 



(7) 
(8) 

(9) 



We will now compute the mass spectrum of each of the models using the generic expressions for the soft masses 
introduced in this section. We will only indicate the new contributions arising from the interaction terms in addition to the 
usual contributions given in Eq. [3] For each model we will indicate the R-parity and the flavor charges for the messenger 
sector. The corresponding charges for the MSSM multiplets are given in Table Q] We will assume that wherever the 
interaction terms include two multiplets, a messenger parity is imposed under which the messenger sector is odd while 
the multiplets of the visible sector are even. 



3.1 Only Singlets 

In these scenarios one has to assume that the messenger sector in addition to the singlets, also has the usual messenger 
fields charged under the MSSM gauge group. As a definite choice we will assume that along with the singlet the messenger 
sector consists of a single 5 © 5 vector pair of chiral messengers (<f> m <f> m ). However a messenger parity prevents them 
from directly coupling with the visible sector. These spectator messengers will contribute to the soft masses through usual 
gauge interactions according to Eq. [3] 



3.1.1 Model 1 

^H^H^-m '■ Considering that the singlet (l m ) is even under R-parity, one obtains the following messenger-matter mixing 
superpotential term iPTD . 

W mt = XH u H d S m . (10) 
The new contributions to the soft masses can be read off from Eqs.[5]and[8] They are given by, 
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(11) 



where a\ i — Xf/Air and the subscripts have their usual meaning. We will follow this convention through out the paper. 
Unfortunately as commented in ifTTl . one can anticipate an anomalous contribution to the (i—B^ parameters in this model. 
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3.1.2 Model 2 



5H51 m : The other possibility here is to assume that the messenger field is odd under the R-parity. However in this case 
we get contributions to the neutrino mass arising through type I see-saw mechanism (26). This puts a lower bound on the 
messenger scale at ~ 10 10 GeV. A way to evade this to is consider another singlet field S and impose non-trivial flavor 
charges: 



U(l) F (l m , l m ,5 m ,5 m ,X) = (-p,-a + p,-(a + q),q,a). 



(12) 



We obtain the following superpotential through which the singlet messenger field interact with the MSSM multiplets, 

W int = XH u L 3 S m . (13) 

We made a field redefinition so that only the third family can interact with the messenger field. We considered a and q 
to be large positive numbers so that all the non-renormalizable terms in the messenger-matter mixing sector are highly 
suppressed. As is clear from the above expression, the field S m must carry a unit lepton number. The new contributions 
to the soft masses from the given interaction term is given by, 
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3.2 Only 5 © 5 

In this section we look at models where the messengers are in the 5 © 5 representations. Depending on the choice of 
symmetries various invariants can be constructed. We will now study the possible terms in turn and compute the new 
contributions to the scalar masses. Again the usual contributions are given by Eq.[3]where we set the Dynkin index d = 1. 
Unless mentioned we will assume that the number of generations of messenger is one. 



3.2.1 Model 3 



105H5 m : Let us consider that the messengers are odd under the R-parity and the following flavor charges are imposed, 

U(l) F {5 m ,5 m ,X) = (-a,0,a). (15) 



We obtain the following interaction term in the superpotential, 



Wint = XQQ'iD'^Hd + XeE^H^Hd- 



(16) 



The absence of messenger parity allows the operator 5h 5 m X to be consistent with all other symmetries imposed. This 
can be absorbed in the superpotential above by a basis change. However there is still a loop level mixing between 5h 
and 5 m . This can be rotated away at the one loop order but contributes non trivially to two loop corrections of the soft 
spectrum. The contributions of these mixing terms to the trilinear coupling are numerically negligible compared to the 
unsuppressed one loop contributions and thus can be neglected. We extend the analysis in |[T2| by including all the two 
loop contributions to the scalar masses and considering the effect of one loop mixing between the messenger and matter 
multiplets. The expressions are given by, 
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3.2.2 Model 4 

10105 m : If we consider the messengers (5 m and 5 TO ) are even under the R-parity we can have two different invariants 
10 10 5 m and 10 5 5 m that are possible. However depending on the assignment of the flavor charges, one or the other 
might may become more dominant. We will consider by turn the two extreme scenarios where only one of the invariants 
dominates. Considering the flavor charges, 



U(l) F (5 m ,5 m ,X) = (0,o, -o), 
the dominating part of the superpotential is given by, 

W mt = \ q Q 3 U^H™ + \ u UiD c m El 
The contributions to the soft scalar masses are given by, 
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The lack of messenger parity in this case can lead to tree level mass terms of the form m'5#5 m . It is expected that 
m' ~ fx, where /i is the usual dimensionful parameter in the MSSM superpotential. This mass cannot be suppressed 
without suppressing the operators in Eq. [19] The phenomenology of this model though not identical, closely resembles 
the scenario studied in |[T31 . 

3.2.3 Model 5 

1055 m : This is the other possibility with even R-parity for the messengers. One can distinguish it from Model 4 by 
considering a different flavor symmetry given by, 



U(l) F (5 m ,5 m ,X) = (a + p, -p,-a). 



(21) 



7 



However note that flavor charges assigned in Table[T]do not distinguish between 5 2 and S3. So this is not enough to ensure 
that only third generation will couple strongly with the messenger sector. This leads to considerable contributions to 
FCNC that constraints the messenger scale. An additional complication is related to the mass terms of the form m'5jj5 m 
where m! ~ e p /i. Worse, the interaction can lead to rapid proton decay and suppression of the order of e 5 is not strong 
enough to comply with present experimental bounds. However one can consider flavor symmetries which are uncorrected 
to the origin of the fermion mass hierarchy that can enable the required suppression. We adopt the paradigm that this is 
possible, assuming this we can write down the superpotential as, 



Wim = \ q Q 3 D c 3 HX + A1Q3L3D: 
The contributions to the soft masses can be written as, 
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(23) 



3.2.4 Model 6 



105 m 5 m : We expand the R-symmetry from Z 2 to Z\. The parity of the MSSM multiplets are still given by Table [T] 
whereas the messenger sector now has the following parity, 



Rp{5m> 5 m ) — (^5 m , 25m). 

On top of this we impose the following flavor charges, 

U(l) F (5 m ,5 m ,X) = (-a,0,a), 
This ensures that we have the following dominant superpotential, 

W mt = \QzD c m H™. 

The new contributions are given by, 
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A t = A b = --^A. 

47T 



3.2.5 Models 7 & 8 

105j m 5km : Sticking to the same invariant as Model 6 an interesting scenario develops when we expand the number of 
messenger generations to two. The drastic change is more than a simple duplication of the results in Model 6. The form 
of the Lagrangian and thus the soft spectrum, depends on the symmetries we impose on the theory. In this regard we will 
discuss two slight variants: 



• Model 7: We consider a Z 4 R-parity with, 
and flavor charges, 

U(l) F (5j m ,5 jm ,X) = (-a,0,a). 

• Model 8: We consider the following Z 2 R-parity, 

Rp{^>lm: 5lrm 52m, 52m) = ( 5lmi 5i m , 52m 5 52m)> 

and the following flavor symmetry: 

^(l)F(5im,5i m ,5 2 m,52m,^) = (o - b, b, a + b, -b, -a) . 



(28) 
(29) 

(30) 

(31) 



With these symmetries we obtain the following form of the superpotential for Model 7, 

W int = (\ 11 Q 3 Dl m H 1 / n +\ 22 Q 3 D c 2m H* m )+\ 12 Q 3 Dt^ 

(32) 

The new contributions to the soft spectrum for Model 7 are given by, 
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The corresponding superpotential and the soft spectrum for Model 8 can be obtained by setting An = A 22 = 0. 



3.3 Only 10 © 10 

In this section we will collect the possible messenger-matter interaction terms possible assuming that the messengers are 
in a vector like representation of 10 © 10. 

3.3.1 Model 9 

1010 m 5H + 10 m 10 m 5H : Consider that the messenger fields are odd under R-parity (Rp(lO m ) = — 10 m ) and have the 
following flavor charges, 

l7(l) a (10 m ,10 m ,X) = (0,-o, a) (34) 

The invariant 10 m §2,3 5h has a coupling which is suppressed by the flavor factor and is at least as small as the A r 
Yukawa. We thus neglect it from the discussion and obtain the following interaction superpotential lfl2l . 



Wi, 



\QzU, n H u + \ u Q m U^H u + \hQ m U^ n H u . 



(35) 



The lack of messenger parity in this scenario again manifests into one loop mixing between messenger and matter multi- 
plets. The one and two loop contributions to the soft masses including the mixing effects are given by, 
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(36) 



3.3.2 Model 10 

10 m 10 m 5h : The other alternative is to consider that the messengers are even under R-parity (i?p(10 m ) = 10 m ). In 
this case the messengers can only couple to the Higgs multiplets in the MSSM sector and thus we select the flavor charges 
of the messenger sector to be zero. The corresponding superpotential is given by IfOl , 



Wi n t = XQmU^Hu, 

and the contributions to the soft masses are given by, 
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3.4 Mixed messenger models 



It is possible that the messenger sector is composed of messenger fields that are in different complete representations of 
the GUT group SU(5). Actually this is implicitly assumed in Models 1 and 2. In that case we can have scenarios where 
more than one of them simultaneously interact with the visible sector. There are a large number of possibilities in this 
class, mainly coming from a combination of two or more models already discussed, possibly augmented by some new 
terms. A study of all these models are beyond the scope of this study and we will restrict ourselves to models that lead to 
entirely new messenger-matter interaction terms. 



3.4.1 Model 11 



Let us consider a scenario where messengers in 1 and 5 representations interact simultaneously with the MSSM multiplets. 
The only new interaction term that can be envisaged in this case is 5h 5 m l m . In order to prevent other terms from 
showing up we impose the following symmetries. Again we conjecture the existence of a second singlet S m . We impose 
the following assignment of R-parity, 



Rpi^mi 5m, lm, lm) — ( 5m, 5m, lm? l-m)i 

and the following flavor charges, 

U(l) F (5 m ,5 m ,l m ,l m ,X) = (-a-b,b,-b,-a + b,a) with a > b. 
We obtain the following unsuppressed terms in the superpotential |fl~3), 

Wi n t = HuH^Sm, 

and the contributions to the soft masses are given by, 
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3.4.2 Model 12 



The messenger sector can very well be made up of singlets and 10 ® 10 multiplets. If both are allowed to couple we can 
have the following new mixed invariant 1010 m l m . In order for this to be possible we can assume that the messenger 
fields have the following R-parity, 



_Z?J)(10 m , 10 m , lm, lm) ( 10 m , 10 m , I771; lm)) 

and the following flavor charges, 

U(l) F (10 m , 10 m , l m , lm, X) = (a-b,b,-b,a + b,-a). 
The leading terms in the superpotential are given by, 

Wi n t — \Q3Qr11Srn + ^uU^U C S m + XeE^E^Sm, 

the new contributions to the soft scalar masses are given by, 
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3.4.3 Models 13 & 14 



Consider scenarios where the messenger sector is made up of vector pairs of 5© 5 and 10© 10. There are two distinct new 
interaction terms that can arise other than combinations of the models already studied earlier. We are going to consider 
these two model one by one, 

1010 m 5 m : One can motivate this by considering the following symmetry arrangements. Consider the following 
R-parity, 

i?p(5 m ,5 m ,10 m ,10 m ) = (-5 m ,-5 m ,10,T0 m ), (48) 

and the following flavor charges, 

J7(l) F (5 TO ,5 m! 10 m ,T0 m ,X) = (6, a -6, -b, a + b, -a), (49) 
with a >> 1 as always. Thus the part of the superpotential that remains unsuppressed by the flavor factors is given by, 

W mt = \Q 3 U 7 c n H™, (50) 
and the new contributions in this case are given by, 
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10 m 5 m 5H : In this case consider the following R-parity, 

Rp{^>mi 5mj 10, 10 m ) = (5 m , 5 m , 10 m ,10 m ), 

and the following flavor charges, 

U(l)F(5 m , 5 m , 10 m , W m ,X) = (b,a-b,a + b, -6, -a). 
The corresponding superpotential is given by, 

Wi n t = XqQ m D^ n H u + A e i? TO i7™-ff u , 
which results in the following new contributions to the soft masses, 
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4 Comparison and Results 

In this section we will compare the models introduced earlier. The ability of these models to raise the Higgs mass, without 
requiring a large stop mass, is through a sizable top trilinear coupling A t . However it is clear from the expressions for the 
soft spectrum given in Eqs. [5] and [8] that the terms that lead to a large A t are also responsible for tachyonic contribution 
to the soft masses. The size of the trilinear coupling is thus constrained by the condition that all sfermion masses should 
remain positive and that proper radiative electroweak symmetry breaking (REWSB) should take place. In Table|2]we list 
the models along with the sfermions that are most vulnerable to the tachyonic contributions. Assuming that the one loop 
tachyonic contributions are negligible owing to the extra suppression from the messenger breaking scale, it is possible 
to put upper limit on the size of the trilinear couplings. We exhibit these upper limits for models wherever they are 
analytically possible by relating them to the gluino mass, which can be considered as the representative of the soft masses. 
Some observations become apparent from the table, 

1 . The right handed sleptons get the smallest contribution from gauge interactions. Thus in models where they receive 
negative contributions from the messenger-matter interactions they become susceptible to turn tachyonic. This 
constrains the size of the trilinear couplings in these models and thus effects their efficiency in raising the Higgs 
masses. 

2. Radiative electroweak symmetry breaking either requires a large stop mass that can turn one of the neutral Higgs 
eigenvalue negative radiatively or a contribution from the messenger-matter interactions that drives it to a negative 
value. An analytic study of each of these models in terms of electroweak symmetry breaking is difficult and we will 
rely on a numerical simulation of these models for this purpose. 

The discussion above is confined to the spectrum at the supersymmetry breaking scale which already gives us an insight 
to models in terms of their ability to generate large trilinear couplings and thus ease the fine-tuning in the models. We 
point out that a generation of a large trilinear is not enough to alleviate the problem with the Higgs mass. In principle 
the subsequent renormalization group evolution to low energies can wash out the trilinear coupling generated at the 
supersymmetry breaking scale due to the partial cancellation between the Yukawa and gaugino contributions to the (3 
function ll27l . In fact a large A t at high scale is useful if it has a sign opposite to that of the gaugino masses. In the 
class of models studied here, this is naturally achieved as is evident from Eq. [3] and Eq.|5] The sign difference between 
the trilinear couplings and the gaugino masses can be traced back to the overall negative factor for fermionic loops in 
the leading contribution to the trilinear coupling. This generic feature of this class of models is aided by our choice of a 
low scale of supersymmetry breaking which reduces any potential harmful effect of the renormalization group running by 
shortening the range. 

Next we perform an extensive numerical study of each of the above models to make a numerical comparison between 
them. The numerical procedure that is followed for each of these models is described below: 

1. We consider A = F/M and x = A/AI as the two independent parameters that define all the scales in the theory. 
In order to scan over the parameter space of the models we vary A between 4 x 10 4 — 9 x 10 5 GeV. We consider 
two different values of x ~ .01, .1 to include scenarios where the one loop contributions become insignificant and 
comparable to the corresponding two loop contributions respectively. 
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Table 2: We summarizer the main theoretical limit on the generated top trilinear coupling where they are analytically possible. We 
also list the smallest value of the lightest stop (mi \min) for which the Higgs mass could be boosted to the range 123 — 127 GeV as 
obtained from our numerical scanning. Finally the NLSP for each model is listed. We also list the references for models which have 
been studied previously in the literature. * indicates that these models are slight variants of the references given or we have updated 
the older analysis to include all two loop corrections. Also note that the eventually the 'unconstrained' are limited by renormalization 
group running effects. 
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2. We use the known Standard Model fermion masses at the weak scale and tan j3 as the boundary conditions and use 
one loop renormalization group equations to determine the known coupling constant at the messenger scale defined 
by M mess = A/x. 

3. For given values of A, x, MSSM couplings and a given set of the new messenger-matter Yukawas we generate the 
soft spectrum at the messenger scale. 

4. We then use one loop renormalization group equations 11271 to determine the weak scale spectrum. Within the frame- 
work of pure gauge mediation the /i term in the MSSM superpotential cannot be generated. Note that messenger- 
matter interactions can in principle lead to generation of a /j, term at the messenger scale, see H13II28I . However all 
the models studied in this paper except Model 1 do not give rise to a (i terrr0. In our numerical study we will simply 
impose the matching condition at the weak scale for the potential minima and set the weak scale /i parameter from 
the following relation, 

m| ~ -2(m 2 Hu + |/i| 2 ) + ^2p( m H d - ™ 2 hJ- (56) 

5. Finally we use SuSpect ll29l and micrOMEGAS [30. 31] to determine correct electroweak symmetry breaking, 
diagonalize the sparticle mass matrices and determine weak scale observables. 

6. We scan over all the new messenger-matter couplings A; between .05 — 1.5. For every choice of the Yukawa 
couplings we repeat the above procedure to scan over the parameter space of the models. 

In FigureQ]we correlate the maximum attainable Higgs mass to the lightest stop mass for the various Models discussed 
in Section [3] We take care to sum over all model parameters and extract the outer envelop of the valid model points to 
obtain the displayed curves. The horizontal band corresponds to the presently preferred range of Higgs mass measured 
at the LHC. This allows a direct comparison between the models introduced in the previous section. To reduce clutter 
we have not displayed the plot for Model 8 which is a subset of the Model 7 in the figures. As anticipated from Table [2] 
models where the size of the trilinear coupling is constrained at the messenger scale show marginal improvement over the 
pure GMSB models. Considerations of radiative electroweak symmetry breaking also enter the game. The table indicates 
the smallest stop mass for which we obtain a Higgs in the desired mass range for every model. A discussion about the 
correlation between the constraints on the trilinear couplings and the lightest stop mass is complicated by the effects of the 
renormalzation group running. However a couple of comments regarding some of the apparent contradictions in Table |2] 
are in order: (i) Note that though the limit on the trilinear coupling is more relaxed in Model 4 as compared to Model 3, 
the former recieves large positive two loop contributions to M? and M"~ Jc thus making the stop spectrum heavier. This 

should be contrasted with the negative constribution to in Model 3. Also an aiding factor is the negative contribution 
to Mfj^ in Model 3 that facilitates a easier REWSB. (ii) While the trilinear couplings in Models 10 and 1 1 have the same 
limits they yeild very different low energy spectrums. This is again related to the difference in the contributions to Mjj 
resulting a more stringent constraint from REWSB for Model 1 1 as compared to Model 10. Note that Model 14 effectively 
combines Models 10 and 11, however a larger Dynkin index makes the spectrum relatively heavier. 

Clearly from this Models 5, 7 ,9, and 10 can be identified as most promising in terms of boosting the Higgs mass to 
experimentally favored range without admitting too much fine-tuning. In Figure [2] we show the allowed region for the 
Models 5, 7 ,9 and 10 where the Higgs mass is between 123 — 127 GeV, in the parameter space of the lightest stop and 
the gluino masses. These plots can be directly translated to a measure of naturalness of the models by relating them to the 
Barbieri-Giudice 11321 fine-tuning parameter using the approximate relation [[331 . A ~ O(l)10i/33(AtMs/650GeV) 2 , 
where t = \og[M mess / Mz], At is the top Yukawa coupling and Ms = ^/m ^ m t ~ 2 . A recent study in pure mGMSB 
models 11341 indicates that a stop mass beyond 4 TeV would be required to obtain a 123 GeV Higgs that tantamounts to a 
considerable fine-tuning. We note that the benchmark points given in Table [3] implies an improvement in the fine-tuning 
at A/ A m cMSB — 0.04 — 0.14. Thus an order of magnitude improvement in fine-tuning can be achieved through the 
messenger-matter interactions discussed in this paper. 

Within the framework of GMSB models, the collider constraints on the sparticle masses are rather model sensetive, 
crucially depending on the NLSP [35.1. Interestingly we observe that among the preferred models only one has a nutralino 

3 However generation of the fj, term through gauge mediation can be incorporated in these models through the usual extension to a Zz NMSSM like 
scenario [13]. 
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Table 3: Representative spectrum for the chosen models. All the masses except the gravitino mass are given in the GeV unit. Here 
tan ft = 10 and x = 0.1 The spectrum is chosen conservatively to be consistent with latest experimental bounds from direct collider 
studies and low energy observables. 
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NLSP while the rest have slepton NLSPfl The search strategies and hence the consequent constrains are different for 
the two scenarios. In case of a bino NLSP one expects a 27 + E™ lss signal, the limits were discussed in a fairly model 
independent manner through simplified models in 1371 . The limits presented for the natural SUSY (light stop) scenario 
indicate trig > 1.1 TeV and > 700 GeV. Subsequent updates [ 38 , 39 1 marginally enhance the exclusion on the gluino. 
The models with slepton NLSP are constrained from SS/OS dileptons + E™ lss and r + leptons + jets + E™ l3S 
searches. Conservative limits around m~ g > 800 GeV were suggested in 11351 . Recent updates from the CMS, see for 
example |40*-42] and the ATLAS, see l43l will certainly increase this bound. However we could not find a comprehensive 
model independent study of this scenario in terms of simplified models in the literature. 

In Table [5] we present indicative benchmark points for the four preferred models. Clearly the limits on the neutralino 
NLSP scenario already push the fine-tuning of the models beyond 10. Let us also point out that recent study on the 
gravitino phenomenology indicates that a cosmologically safe upper limit gravitino mass might be as low as 16 eV B31 . 
The gravitino mass is given by the usual relation m 3 / 2 = F '/ (\<p\/3AIpi) where F and are defined in Eq.|2]and Mpi 
is the Planck mass. Comparing with the gaugino masses given in Eq.[3]we see that, 

TO 3/2 M 1 

M g ^ Mpi N mess d' 

where M,N mess and d are the messenger scale, messenger multiplicity and the Dynkin index respectively. It is clear that 
the restrictive limits from the gravitino mass can be accomodated in models with higher messenger multiplicity and/or 
high Dynkin index. It might be expected that the large gravitino mass (e.g. the benchmark point for Model 5 in Table [3) 
for models with low Dynkin index can be handled by considering higher messenger multiplicity. A detailed study of the 
effect of varying the number of messengers is beyond the scope of this paper. 

Note that most of the discussion in this paper corresponds to this light gravitino scenario. However as the gravitino mass 
crosses ~ 50MeV the NLSP becomes increasingly stable leading to stronger constraints from searches for long living 
particles at the LHC. For example a quasi-stable stau has the constraints > 223 GeV on its mass [4411 . Consequently 
this would rule out most of the low fine-tuned regions in our models providing an added motivation to consider low scale 
of supersymmetry breaking. 



5 Conclusion 



Supersymmetric scenarios with viable UV complete supersymmetry breaking sectors are severely challenged by the mea- 
sured Higgs mass at the LHC. Several extension of these models that conjecture the existence of extra exotic particles 
have been proposed in the literature to address this issue [46-511] . Within gauge mediated models one can evade these 
constraints economically, by considering messenger-matter interactions. Some of these models have also been studied in 
the recent past , ifTTTfTBI . In the present paper we study a class of these models that can be embedded in SU(5) GUT group. 
We find R-parity and flavor symmetries provide an organized way to study these models. Interestingly the flavor symme- 
try can be tied to the origin of the observed SM fermion masses and mixing through a Froggatt-Nielsen like framework. 
We construct all the possible invariant terms with the messengers in the 1, 5 and 10 dimensional representations of the 
SU(5) group. Some models or their close variants studied in the literature show up naturally within this framework while 
a whole set of new models are predicted. However it should be emphasized that the choice of flavor symmetries presented 
here is not unique. They are presented mostly in the spirit of proof of principle. It might as well be that the form of the 
interaction superpotential discussed here arises from a completely different underlying flavor structure. Nevertheless the 
phenomenological features studied in the second half of the paper are independent of these assumptions. 

A detailed numerical study of the relevant parameter space of the different scenarios is carried out to compare these 
models. We identify models that can effectively raise the Higgs mass to the favoured range without admitting too much 
fine-tuning. We find many of the generic and specific SUSY searches at LHC imply considerable constraints on these 
models. Incidentally we observe that all but one of the preferred models have a slepton NLSP. However we could not 

4 As a small digression it is curious to note that a light stau can lead to enhanced contribution to the Br(h — » 77) rate through loop contributions in 
the large tan /3 regime 1361 and can explain the slight excess over SM prediction indicated by the present experimental data. 
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Fi gure 1 : The maximum Higgs mass for the different models is plotted against the lightest stop mass. In all the models the new Yukawa 
couplings are varied between 0.1 — 1.5, tan B between 5 — 50, A between 2 x 10 4 — 9 x 10° and two different choices of x = .01 & .1 
were considered. The region below each curve is accessible to the corresponding model. 'GMSB' in (a) represents the mGMSB model 
with a single messenger pair in 5 © 5 of SU(5). 



locate a detailed model independent study of the LHC constraints on this scenario in the literature. Most studies are 
restricted to benchmark points or related to details of the entire spectrum. Considering that many other possible natural 
SUSY models within the context of GMSB could lead to a slepton NLSP, a detailed study of this scenario is highly 
anticipated. 
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Figure 2: Allowed region plot of the preferred models which correspond to nth between 123 — 127 GeV in the plane of the lightest 
stop mass and the gluino mass. 

References 

[1] G. Aad et al. [ATLAS Collaboration], "Observation of a new particle in the search for the Standard Model Higgs 
boson with the ATLAS detector at the LHC," Phys. Lett. B 716 (2012) 1 [arXiv:1207.7214 [hep-ex]]. 

[2] S. Chatrchyan et al. [CMS Collaboration], "Observation of a new boson at a mass of 125 GeV with the CMS 
experiment at the LHC," Phys. Lett. B 716 (2012) 30 [arXiv: 1207.7235 [hep-ex]]. 

[3] P. Draper, P. Meade, M. Reece and D. Shih, "Implications of a 125 GeV Higgs for the MSSM and Low-Scale SUSY 
Breaking," Phys. Rev. D 85 (2012) 095007 [arXiv:l 1 12.3068 [hep-ph]]. 

[4] G. F. Giudice and R. Rattazzi, "Theories with gauge mediated supersymmetry breaking," Phys. Rept. 322, 419 
(1999) [hep-ph/9801271]. 

[5] K. Intriligator and M. Sudano, "General Gauge Mediation with Gauge Messengers," JHEP 1006 (2010) 047 
[arXiv: 100 1.5443 [hep-ph]]. 

[6] B. Bajc, S. Lavignac and T. Mede, "Supersymmetry Breaking Induced by Radiative Corrections," JHEP 1207 (2012) 
185 [arXiv: 1202.2845 [hep-ph]]. 

[7] R. Barbieri, L. J. Hall and A. Strumia, "Violations of lepton flavor and CP in supersymmetric unified theories," Nucl. 
Phys. B 445 (1995) 219 [hep-ph/9501334], 

[8] R. Barbieri, L. J. Hall and A. Strumia, "Hadronic flavor and CP violating signals of superunification," Nucl. Phys. B 
449 (1995) 437 [hep-ph/9504373], 

[9] N. Arkani-Hamed, H. -C. Cheng and L. J. Hall, "Flavor mixing signals for realistic supersymmetric unification," 
Phys. Rev. D 53 (1996) 413 [hep-ph/9508288]. 

[10] M. Y. .Khlopov and A. D. Linde, "Is It Easy to Save the Gravitino?," Phys. Lett. B 138 (1984) 265. 



19 



[11] Z. Kang, T. Li, T. Liu, C. Tong and J. M. Yang, "A Heavy SM-like Higgs and a Light Stop from Yukawa-Deflected 
Gauge Mediation," Phys. Rev. D 86, 095020 (2012) [arXiv: 1203.2336 [hep-ph]]. 

[12] A. Albaid and K. S. Babu, "Higgs boson of mass 125 GeV in GMSB models with messenger-matter mixing," 
arXiv:1207.1014 [hep-ph]. 

[13] N. Craig, S. Knapen, D. Shih and Y. Zhao, "A Complete Model of Low-Scale Gauge Mediation," arXiv: 1206.4086 
[hep-ph]. 

[14] M. Abdullah, I. Galon, Y. Shadmi and Y. Shirman, "Flavored Gauge Mediation, A Heavy Higgs, and Supersymmetric 
Alignment," arXiv: 1209.4904 [hep-ph]. 

[15] J. L. Evans, M. Ibe and T. T. Yanagida, "Relatively Heavy Higgs Boson in More Generic Gauge Mediation," Phys. 
Lett. B 705, 342 (2011) [arXiv: 1107.3006 [hep-ph]]. 

[16] S. P. Martin and J. D. Wells, "Implications of gauge-mediated supersymmetry breaking with vector-like quarks and 
a 125 GeV Higgs boson," Phys. Rev. D 86, 035017 (2012) [arXiv: 1206.2956 [hep-ph]]. 

[17] L. Calibbi, L. Ferretti, A. Romanino and R. Ziegler, "Gauge coupling unification, the GUT scale, and magic fields," 
Phys. Lett. B 672, 152 (2009) [arXiv:0812.0342 [hep-ph]]. 

[18] P. Byakti and D. Ghosh, "Magic Messengers in Gauge Mediation and signal for 125 GeV boosted Higgs boson," 
Phys. Rev. D 86, 095027 (2012) [arXiv: 1204.0415 [hep-ph]]. 

[19] C. D. Froggatt and H. B. Nielsen, "Hierarchy of Quark Masses, Cabibbo Angles and CP Violation," Nucl. Phys. B 
147(1979) 277. 

[20] C. A. Savoy and M. Thormeier, "Exotic particles below the TeV from low scale flavour theories," JHEP 1103 (201 1) 
128 [arXiv: 1003.2090 [hep-ph]], 

[21] O. J. P. Eboli, C. A. Savoy and R. Zukanovich Funchal, "A Rationale for Long-lived Quarks and Leptons at the LHC: 
Low Energy Flavour Theory," JHEP 1202 (2012) 123 [arXiv:1112.5108 [hep-ph]], 

[22] T. S. Ray, H. de Sandes and C. A. Savoy, "Gluino, Wino and Higgsino-Like Particles without Supersymmetry," Phys. 
Lett. B 712 (2012)401 [arXiv:l 1 12.6180 [hep-ph]]. 

[23] S. Dimopoulos and G. F. Giudice, "Multimessenger theories of gauge mediated supersymmetry breaking," Phys. 
Lett. B 393 (1997) 72 [hep-ph/9609344]. 

[24] G. F. Giudice and R. Rattazzi, "Extracting supersymmetry breaking effects from wave function renormalization," 
Nucl. Phys. B 511, 25 (1998) [hep-ph/9706540] . 

[25] Z. Chacko and E. Ponton, "Yukawa deflected gauge mediation," Phys. Rev. D 66, 095004 (2002) [hep-ph/01 12190]. 

[26] R. N. Mohapatra and P. B. Pal, "Massive neutrinos in physics and astrophysics. Second edition," World Sci. Lect. 
Notes Phys. 60, 1 (1998) [World Sci. Lect. Notes Phys. 72, 1 (2004)]. 

[27] S. P. Martin, "A Supersymmetry primer," In *Kane, GL. (ed.): Perspectives on supersymmetry II* 1-153 [hep- 
ph/9709356]. 

[28] A. Delgado, G. F. Giudice and P. Slavich, "Dynamical mu term in gauge mediation," Phys. Lett. B 653 (2007) 424 
[arXiv:0706.3873 [hep-ph]]. 

[29] A. Djouadi, J. -L. Kneur and G. Moultaka, "SuSpect: A Fortran code for the supersymmetric and Higgs particle 
spectrum in the MSSM," Comput. Phys. Commun. 176 (2007) 426 [hep-ph/021 1331]. 

[30] G. Belanger, F. Boudjema, P. Brun, A. Pukhov, S. Rosier-Lees, P. Salati and A. Semenov, "Indirect search for dark 
matter with micrOMEGAs2.4," Comput. Phys. Commun. 182 (201 1) 842 [arXiv: 1004.1092 [hep-ph]]. 

[31] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, "micrOMEGAs: Version 1.3," Comput. Phys. Commun. 
174 (2006) 577 [hep-ph/0405253]. 



20 



[32] R. Barbieri and G. F. Giudice, "Upper Bounds on Supersymmetric Particle Masses," Nucl. Phys. B 306 (1988) 63. 
[33] D. Wright, "Naturally nonminimal supersymmetry," hep -ph/980 1449. 

[34] M. A. Ajaib, I. Gogoladze, F. Nasir and Q. Shan, "Revisiting mGMSB in Light of a 125 GeV Higgs," Phys. Lett. B 
713 (2012) 462 [arXiv: 1204.2856 [hep-ph]]. 

[35] Y. Kats, P. Meade, M. Reece and D. Shih, "The Status of GMSB After 1/fb at the LHC," JHEP 1202 (2012) 115 
[arXiv:l 1 10.6444 [hep-ph]]. 

[36] M. Carena, S. Gori, N. R. Shah, C. E. M. Wagner and L. -T. Wang, "Light Stau Phenomenology and the Higgs 77 
Rate," JHEP 1207 (2012) 175 [arXiv: 1205.5842 [hep-ph]]. 

E3 

[37] J. Barnard, B. Farmer, T. Gherghetta and M. White, "Natural gauge mediation with a bino NLSP at the LHC," 
arXiv: 1208.6062 [hep-ph]. 

[38] G. Aad et al. [ATLAS Collaboration], "Search for diphoton events with large missing transverse momentum in 7 
TeV proton-proton collision data with the ATLAS detector," arXiv: 1209.0753 [hep-ex], 

[39] S. Chatrchyan et al. [CMS Collaboration], "Search for new physics in events with photons, jets, and missing trans- 
verse energy in pp collisions at sqrt(s) = 7 TeV," arXiv:121 1 .4784 [hep-ex]. 

[40] S. Chatrchyan et al. [CMS Collaboration], "Search for new physics in events with same-sign dileptons and 6-tagged 
jets in pp collisions at = 7 TeV," JHEP 1208 (2012) 1 10 [arXiv: 1205.3933 [hep-ex]], 

[41] S. Chatrchyan et al. [CMS Collaboration], "Search for new physics with same-sign isolated dilepton events with jets 
and missing transverse energy," Phys. Rev. Lett. 109 (2012) 071803 [arXiv:1205.6615 [hep-ex]], 

[42] S. Chatrchyan et al. [CMS Collaboration], "Search for new physics in events with opposite-sign leptons, jets, and 
missing transverse energy in pp collisions at sqrt(s) = 7 TeV," arXiv: 1206.3949 [hep-ex]. 

[43] G. Aad et al. [ATLAS Collaboration], "Search for Supersymmetry in Events with Large Missing Transverse Mo- 
mentum, Jets, and at Least One Tau Lepton in 7 TeV Proton-Proton Collision Data with the ATLAS Detector," 
arXiv:1210.1314 [hep-ex]. 

[44] S. Chatrchyan et al. [CMS Collaboration], "Search for heavy long-lived charged particles in pp collisions at sqrt(s)=7 
TeV," Phys. Lett. B 713 (2012) 408 [arXiv: 1205.0272 [hep-ex]]. 

[45] M. Viel, J. Lesgourgues, M. G. Haehnelt, S. Matarrese and A. Riotto, "Constraining warm dark matter candidates 
including sterile neutrinos and light gravitinos with WMAP and the Lyman-alpha forest," Phys. Rev. D 71 (2005) 
063534 [astro-ph/0501562]. 

[46] T. Moroi, R. Sato and T. T. Yanagida, "Extra Matters Decree the Relatively Heavy Higgs of Mass about 125 GeV in 
the Supersymmetric Model," Phys. Lett. B 709 (2012)218 [arXiv:1112.3142 [hep-ph]]. 

[47] U. Ellwanger and C. Hugonie, "Higgs bosons near 125 GeV in the NMSSM with constraints at the GUT scale," 
Adv. High Energy Phys. 2012 (2012) 625389 [arXiv: 1203.5048 [hep-ph]]. 

[48] K. J. Bae, T. H. Jung and H. D. Kim, "125 GeV Higgs as a pseudo-Goldstone boson in supersymmetry with vector- 
like matters," arXiv: 1208.3748 [hep-ph]. 

[49] B. Kyae and J. -C. Park, "A Singlet-Extension of the MSSM for 125 GeV Higgs with the Least Tuning," 
arXiv: 1207.3 126 [hep-ph]. 

[50] G. Bhattacharyya and T. S. Ray, "Pushing the SUSY Higgs mass towards 125 GeV with a color adjoint," 
arXiv:1210.0594 [hep-ph], 

[51] M. Endo, K. Hamaguchi, K. Ishikawa, S. Iwamoto and N. Yokozaki, "Gauge Mediation Models with Vectorlike 
Matters at the LHC," arXiv:1212.3935 [hep-ph]. 



21 



